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On  the  Pressure  Drop  Through  a  Uniform  Snow  Layer 


YIN-CHAO  YEN 


INTRODUCTION 

Fluid  flow  through  a  bed  of  contiguous  particles 
is  distinct  from  the  flow  through  a  conduit.  Since 
its  passage  is  between  the  particles  of  the  bed,  the 
flow  is  dependent  upon  the  porosity  of  the  bed  or 
column  and  diameter,  sphericity,  orientation  and 
roughness  of  the  particles.  The  actual  linear  veloc¬ 
ity  of  the  fluid  through  the  passages  in  the  porous 
medium  may  be  expressed  in  terms  of  “superfi¬ 
cial”  velocity  (computed  as  the  rate  of  flow  of 
fluid  through  the  entire  undisturbed  cross-sectional 
area  of  the  bed). 

The  most  frequently  used  model  for  volume- 
averaged  flow  through  a  porous  medium  is 
Darcy’s  flow  model.  In  one-dimensional  flow,  it 
can  be  expressed  by 


vrSN 

SA 


(3) 


and 


K  = 


b> 

\2(a  +  b) 


(4) 


where  N  is  the  number  of  tubes  contained  on  a 
cross  section  of  area  A. 

Kozeny  (1927),  on  the  other  hand,  considered 
the  porous  medium  as  a  collection  of  solid  spheres 
and  arrived  at  the  following  permeability  expres¬ 
sion: 


K  - 


d2t' 

(l-0! 


(5) 


where  u  =  velocity  component  in  .v-direction 
K  =  permeability  of  the  medium 
n  =  viscosity 
p  =  pressure. 

For  three-dimensional  flow  in  the  presence  of  a 
gravitational  acceleration  vector  g  (g„,  gy,  gz),  eq  1 
can  be  written  as 

v  =  —  (-  Vp  +  eg)  (2) 

n 

where  v  is  the  velocity  vector,  p  is  density  and  K, 
usually  expressed  in  units  of  length  squared  (i.e. 
m\  cm1,  or  ft1),  can  be  determined  by  measuring 
the  pressure  drop  and  the  flow  rate.  Considering 
the  pores  as  a  bundle  of  parallel  capillary  tubes  of 
radius  r„  or  the  pores  as  a  stack  of  parallel  capil¬ 
lary  fissures  of  width  b  and  fissure-to-fissure  spac¬ 
ing  a  +  b.  Bear  (1972)  derived  the  following  ex¬ 
pressions: 


where  <  is  the  porosity  (the  ratio  of  void  volume  to 
the  total  volume  of  the  porous  medium),  and  d  is 
the  diameter  of  the  spheres.  Darcy’s  flow 
model  is  believed  to  be  valid  as  long  as  the  so-called 
local  pore  Reynolds  number  Re,  based  on  local 
volume-averaged  speed  [i.e.  (w1  +  v1  +  tv1)'1]  and 
K 1/1  are  smaller  than  one.  For  Re,  greater  than  one, 
modification  of  this  model  has  been  proposed  by 
Fcrschheimer  (1901)  as 

=  (6) 

For  three  dimensions  and  in  the  presence  of  body 
acceleration,  this  modification  becomes 

BeK  K  .  rr  \  n\ 

ft  /< 

Ward  (1964)  suggested  the  value  of  B  as 

B  =  -0.55K'/!.  (8) 

for  Re,  exceeding  approximately  10.  Ergum  (1952), 


< 


after  conducting  extensive  measurements  on  gas 
flow  through  columns  of  packed  spheres,  sands 
and  pulverized  coal,  provided  the  correlations  for 
K  and  B  as 


150(1  -f)J 


and 


B 


1.75(1  - t) 

dt' 


(9) 


(10) 


The  loss  due  to  friction  accompanying  the  flow 
through  porous  media  can  be  derived  from  expres¬ 
sions  used  for  fluid  flow  through  conduits  such  as 


2  gcD/wf 

f  L  v!  Hlo/W  /I 


(ID 


where  is  the  mass  fraction  of  a  given  particle 
size  and  D\  is  the  diameter  of  particles  in  each  size 
fraction  taken  as  the  arithmetic  average  of  the 
screen  openings  passing  and  retaining  the  parti¬ 
cles.  The  value  of  FRe  and  Ff  can  be  evaluated  from 
the  graph  given  in  Brown  et  al.  (1953)  and  present¬ 
ed  as  functions  of  t  and  v.  The  porosity  is  closely 
related  to  the  sphericity.  However,  different  por¬ 
osities  are  expected  with  particles  of  the  same 
shape  through  variations  in  spatial  arrangement. 
Consequently  both  porosity  and  sphericity  are  re¬ 
quired  to  define  the  porous  medium.  Although  the 
sphericity  could  theoretically  be  calculated  from 
the  dimensions  of  the  particle,  this  is  often  diffi¬ 
cult  or  virtually  impossible  for  randomly  packed 
beds  of  uniform-sized  particles.  The  sphericity  is 
found  to  be  a  unique  function  of  porosity  with  a 
slight  variation  due  to  the  nature  of  packing  (i.e. 
loose,  normal  or  dense). 


where  / 

8c 

D 

L 

/w, 

u 


frictional  loss  factor 

conversion  factor  (32.17  poundals  per 

pound  in  English  units) 

diameter  of  the  conduit 

length  of  the  conduit 

frictional  energy  loss  per  unit  mass 

surface  roughness 

superficial  velocity,  i.e.  a  linear  veloc¬ 
ity  computed  on  the  basis  of  the  total 
crosssectional  area. 


In  the  case  of  fluid  flow  through  a  porous  medi¬ 
um,  the  Reynolds  number  will  be  modified  to 


Rep  = 


/< 


(12) 


where  Dp  is  the  diameter  of  the  particle,  is  a 
factor  dependent  on  the  value  of  porosity  <  and 
sphericity  v!  (which  is  defined  as  the  area  of  the 
sphere  having  the  same  volume  as  the  particle  di¬ 
vided  by  the  area  of  the  particle)  of  the  medium. 
The  friction  factor /p  is  adjusted  to  include  a  fac¬ 
tor  Ff  also  dependent  on  the  values  of  <  and  v  as 


2gc  Dplwf 
Lv2  Ff 


2f-cDv(-  Ap)fQ 
FfLG2 


(13) 


When  ail  particles  are  of  the  same  size,  the  screen 
size  (Dav g)  may  be  used.  For  mixed  sizes,  Dp  is  the 
mean  surface  diameter  given  by 


EXPERIMENTAL  SETUP 
AND  PROCEDURE 

The  test  apparatus  is  identical  to  the  one  used  in 
determining  the  effective  thermal  conductivity  and 
water  vapor  diffusivity  of  the  snow  (Yen  1962, 
1963).  Figure  1  is  a  schematic  of  the  experimental 
setup.  The  entire  apparatus,  with  the  exception  of 
the  manometer,  the  wet-test  meter  (to  measure  the 
air  flow  rate)  and  the  Leeds  and  Northrup  potenti¬ 
ometer  (placed  in  the  adjacent  corridor  maintained 
at  ordinary  room  temperature),  was  situated  in  a 
refrigerated  room  kept  at  approximately  -20°C. 
Coldroom  air  was  compressed  and  then  passed 
through  a  tank  to  minimize  the  fluctuations  of  the 
air  pressure.  A  pancake-type  pressure  regulator 
was  used  to  obtain  the  desired  flow  rate.  Since  the 
coldroom  could  only  be  maintained  to  within 
±  1  °C,  the  air  was  passed  through  a  constant  tem¬ 
perature  bath  where  the  temperature  could  be 
maintained  at  slightly  above  room  temperature  to 
an  accuracy  of  ±0.05°C. 

The  Dewar  flask,  6.5cm  I.D.,  15.24cm  in  height 
and  made  of  double-walled  glass  tubing  with  the 
annular  space  being  evacuated  and  the  internal 
surface  plated  with  silver,  contained  the  snow 
sample  and  was  enclosed  in  an  aluminum  casing 
equipped  with  a  flange  at  both  ends.  The  casing 
was  made  so  that  the  flask  fitted  tightly  into  it, 
and  the  flask  was  coated  with  a  low-temperature 
silicone  grease  before  assembly  to  prevent  any  pos¬ 
sible  air  leakage.  As  soon  as  a  steady  state  was  es¬ 
tablished,  the  water  temperature  at  the  wet-test 
meter  was  recorded  with  a  Rubican  potentiometer 
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Air  Inlet 

To  Micromanometer  ■* 


Dewar  Bask 


Wood  Bushing 


To  Micromanometer 


♦To  Potentiometer 
Rubber  Gasket 


ubber  Gasket 


Exit  Air  to  Wet¬ 
test  Meter 


Thermocouple  Wire 
to  Potentiometer 


Figure  I.  Schematic  of  the  experimental  setup. 


accurate  to  ±0.05°C  along  with  the  absolute 
barometric  pressure.  The  wet-test  meter  was 
checked  by  the  standard  displacement  method  and 
was  found  to  be  accurate  within  20  cm!  out  of 
3000  cm’  (±0.7%).  The  displacement  capacity  of 
the  meter  is  300  cm’  per  revolution  with  300  divi¬ 
sions  on  the  dial.  The  snow  samples  for  the  experi¬ 
ments  covered  a  range  of  sieve  numbers,  5-10,  JO- 
12,  14-16,  18-20  and  25-30,  with  corresponding 
nominal  snow  particle  diameters  of  0.219,  0.182, 
0.129,  0.092  and  0.065  cm,  respectively.  The  snow 
sample  density  was  determined  by  recording  the 


amount  of  snow  poured  into  the  flask  to  be  flush 
with  the  sample  container  (i.e.  ps  =  weight  of 
snow/volume  of  the  sample  container).  Several 
determinations  were  made  for  all  the  specific  sieve 
fractions  and  an  arithmetic  average  snow  density 
was  determined. 


EXPERIMENTAL  RESULTS 

A  great  number  of  experimental  runs  were  con¬ 
ducted  for  each  snow  density  covered  in  the  study. 


3 


Table  1.  Summary  of  derived  parameters  as  a  function  of  snow  density. 


Snow 
density,  u, 
tg/cm') 

Snow 
porosity,  i 

Snow 
particle 
diameter,  Dp 
< cm ) 

Sphericity 
of  snow 
particle,  i 

Reynolds 
number 
factor,  fKe 

Friction 

coefficient. 

0.377 

0.589 

0.219 

0.55 

43.5 

830 

0.387 

0.578 

0.182 

0.55 

45.0 

950 

0.427 

0.534 

0.129 

0.60 

47.0 

1150 

0.447 

0.512 

0.092 

0.65 

47.0 

1150 

0.472 

0.485 

0.065 

0.70 

47.0 

1150 

Table  1  summarizes  the  derived  parameters  as 
function  of  snow  density.  The  porosity  e  is  evalu¬ 
ated  as 


where  gj  is  the  ice  density,  usually  taking  a  value 
of  0.917  g/cm*.  The  snow  particle  nominal  diam¬ 
eter  Dp  was  evaluated  by  arithmetically  averaging 
the  two  sieve  openings.  The  sphericity  Reynolds 
number  factor  F*  and  coefficient  Ff  are  taken 
from  the  graphs  given  in  Brown  et  al.  (1953).  Since 
the  snow  density  has  only  slightly  more  than  25% 
variation,  there  is  merely  a  weak  variation  at  the 
most  in  values  of  i  and  F*'.  However,  the  values 
of  Ff  were  rather  sensitive  to  the  changes  of  f  and 


\p.  Figures  2a  and  b  show  the  variation  of  pressure 
drop  ( -  Ap) f  vs  the  mass  flow  rate  G  (g/cm!-s).  It 
can  be  seen  that  under  the  combination  of  experi¬ 
mental  conditions  (i.e.  for  the  range  of  snow  den¬ 
sity  and  the  variation  of  the  mass  flow  rate  cov¬ 
ered  in  this  study)  a  linear  relationship  exists  be¬ 
tween  ( -  Ap)f  and  G.  The  coefficients  a  and  b  of 
the  linear  expression 

(-Ap)f  =  a  +  bGx  10-  (16) 

were  determined  from  least-squares  analysis  and 
are  summarized  in  Table  2  along  with  the  correla¬ 
tion  coefficients.  Theoretically,  the  values  of  a 
should  be  zero,  and  small  values  of  a  are  indicative 
of  the  accuracy  of  the  experimental  measurement. 
However,  for  practical  purposes,  one  can  hardly 


a.  For  Qs  =  0.377,  0.387,  0.427  and  0.447  g/cm\ 


Figure  2.  Pressure  drop  due  to  friction  losses  as  a  function  of  mass  flow  rate 
G. 
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b.  For  qs  =  0.472  g/cm\ 
Figure  2  (coni ’d). 


Table  2.  Coefficients  a  and  b  in  (-A p),  = 
a  +  6G  x  lO  "  and  coerrelation  coefficient  y. 


Snow 
density,  gs 

Coefficient 
a  b 

(lif'/cm)  (gf‘S/g) 

Correlation 

coefficient, 

0.377 

3.671  xIO" 

31.267 

0.9492 

0.387 

1.582x10  1 

39.382 

0.8130 

0.427 

3.137x10  ' 

88.425 

0.8681 

0.447 

8.854x  10  ’ 

176.95 

0.9933 

0.472 

8.478  x10  1 

413.20 

0.9504 

*  gf— gram-force. 


with  a  rather  high  correlation  coefficient  of 
0.9688.  Equation  17  gives  a  value  of/p  approach¬ 
ing  that  of  /p  =  65/Rep  as  Rep  approaches  the 
limit  of  the  laminar  flow  region  (say  Rep  =  400) 
recommended  for  randomly  packed  particles  (me¬ 
tallic  or  nonmetallic  cylinders,  rings,  spheres, 
etc.). 

Another  way  to  represent  the  experimental  re¬ 
sults  is  by  expressing  the  friction  factor /p  and  the 
Reynolds  number  Rep,  respectively,  by 


yP  2(1-0 


(18) 


detect  that  the  lines  in  Figures  2a  and  b  do  not  pass 
through  the  origin. 

Table  3  summarizes  the  experimental  data  and 
the  calculated  values  of  Rep  (defined  as  DpGFRe/ 
n)  and  friction  factor /p  [defined  in  eq  15  as  /p  = 
2gcDp(  -  Ap)fQ  / FfLG2].  An  average  value  of  n 
based  on  the  arithmetic  mean  of  the  inlet  and  out¬ 
let  air  temperature  was  used.  The  mass  flow  rate 
was  determined  on  the  basis  of  the  exit  air  temper¬ 
ature  and  the  prevailing  atmospheric  pressure. 
Figure  3  shows  the  variation  of  fp  versus  Rep.  The 
results  of  a  least-squares  analysis  of  the  results  can 
be  expressed  by 


fp  - 


118 

Rep095 


(17) 


and 


Rep  = 


Rep 

^ReO-O  ‘ 


(19) 


The  least-squares  analysis  of  fp  and  Rep  can  be 
represented  by 


fp  - 


Rep 


Re 


<0.9093 

P 


(20) 


but  only  with  a  moderate  correlation  coefficient  of 
0.8385.  Figure  4  shows  a  plot  of  fp  vs  Rep,  and  it 
can  be  clearly  seen  that  the  value  fp  has  a  weak  de¬ 
pendence  on  the  snow  density  or  porosity. 
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t 


Table  3.  Summary  of  experimental  results  and  calculated  data. 


Ox/S  ‘ 
(g/cm‘  s) 

( gf/cm ’) 

Re 

^D^-JVpVx/FfLGV 

e,  =  0.377  g/cmJ,  0p  = 

0.219  cm 

13.11 

0.048 

75.1 

1.28 

14.02 

0.051 

80.3 

1.19 

9.11 

0.032 

52.2 

1.77 

7.57 

0.027 

43.3 

2.16 

15.09 

0.051 

86.4 

1.03 

16.20 

0.053 

92.8 

0.92 

17.80 

0.062 

102.0 

0.90 

19.52 

0.062 

112.0 

0.74 

24.70 

0.085 

142.1 

0.64 

27.09 

0.088 

155.7 

0.55 

25.09 

0.078 

144.2 

0.57 

22.68 

0.070 

130.4 

0.62 

21.47 

0.068 

123.5 

0.68 

19.98 

0.067 

114.9 

0.77 

5  66 

0.019 

32.4 

2.72 

9.70 

0.033 

55.6 

1.61 

11.63 

0.034 

66.6 

1.15 

14.14 

0.041 

81.0 

0.94 

15.82 

0.052 

90.9 

0.95 

19.77 

0.057 

113.5 

0.67 

21.75 

0.068 

124.8 

0.66 

24.02 

0.075 

137.9 

0.59 

21.91 

0.066 

125.9 

0.63 

10.61 

0.031 

60.9 

1.26 

3.71 

0.012 

21.3 

4.00 

6.84 

0.024 

39.1 

2.35 

16.60 

0.058 

95.0 

0.96 

17.70 

0.055 

101.3 

0.80 

20.93 

0.070 

119.8 

0.73 

22.57 

0.075 

129.1 

0.67 

16.02 

0.051 

91.7 

0.91 

14.70 

0.047 

84.0 

1.00 

13.29 

0.042 

75.9 

1.08 

18.36 

0.054 

105.2 

0.73 

6.23 

0.022 

35.7 

2.60 

7.86 

0.024 

45.1 

1.78 

10.70 

0.029 

61.5 

1.16 

12.34 

0.038 

70.9 

1.14 

15.45 

0.047 

89.0 

0.90 

11.70 

0.035 

67.3 

1.17 

4.66 

0.015 

26.7 

3.17 

13.36 

0.037 

76.9 

0.95 

15.54 

0.038 

89.4 

0.72 

18.20 

0.046 

104.9 

0.64 

19.39 

0.051 

111.8 

0.62 

e,  =  0.387  g/cm\  O p  = 

0.182  cm 

9.70 

0.053 

47.8 

1.88 

10.69 

0.056 

52.7 

1.63 

12.92 

0.070 

63.8 

1.40 

28.83 

0.137 

142  7 

0.55 

29.43 

0.142 

145.6 

0.55 

30.48 

0.144 

150.9 

0.52 

23.49 

0.107 

116.2 

0.65 

13.89 

0.060 

68.8 

1.04 

G  x  10  * 
(g/cm ’  s) 

(gf/cm') 

Re 

K>D¥ReO/n) 

fp 

f2gcOJ-lp)lV/rfLC:J 

16.73 

0.069 

82.7 

0.82 

19.43 

0.075 

95.9 

0.66 

20.94 

0.0% 

103.5 

0,73 

29.20 

0.118 

144.5 

0.46 

25.50 

0.105 

126.2 

0.54 

31.16 

0.123 

154.2 

0.42 

20.92 

0.083 

103.4 

0.63 

25.65 

0.081 

125.9 

0.41 

8.77 

0.051 

43.2 

2.21 

9.67 

0.051 

47.6 

1.81 

10.73 

0.061 

52.9 

1.77 

12.59 

0.068 

62.1 

1.43 

30.93 

0.150 

153.0 

0.52 

29.72 

0.135 

147.1 

0.51 

21.50 

0.091 

106.3 

0.66 

20.26 

0.089 

100.2 

0.72 

19.46 

0.083 

%.3 

0.73 

19.33 

0.081 

95.5 

0.72 

18.33 

0.077 

90.6 

0.76 

16.84 

0.072 

83.2 

0.85 

15.94 

0.066 

78.8 

0.86 

27.80 

0.106 

137.6 

0.46 

22.72 

0.088 

112.3 

0.57 

17.50 

0.102 

86.0 

111 

20.93 

0.123 

86.5 

0.94 

23.27 

0.137 

115.0 

0.84 

25.54 

0.145 

126.3 

0.74 

27.27 

0.157 

134.9 

0.70 

11.75 

0.070 

57.8 

1.69 

14.24 

0.083 

70.2 

1.36 

15.41 

0.085 

71.5 

1.19 

17.52 

0.096 

86.2 

1.04 

18.20 

0.100 

89.8 

1.01 

6.39 

0.044 

31.5 

3.60 

8.00 

0.048 

39.4 

2.50 

10.43 

0.063 

42.8 

1.93 

11.89 

0  065 

48.8 

1.53 

3.39 

0.026 

16.7 

7.56 

13.41 

0.076 

66.1 

1.41 

g.  =  0.427  g/cm‘,  l)f  =  0.129  cm 


3.84 

0.050 

14.1 

6.63 

16.95 

0.190 

62.5 

1.29 

9.48 

0.107 

34.9 

2.33 

4.11 

0.050 

15.1 

5.77 

5.00 

0.059 

18.3 

4.61 

2.48 

0.034 

9.1 

10.82 

18.04 

0.2C1 

66.6 

1.21 

10  16 

0.110 

37.4 

2.08 

24.38 

0.246 

90.1 

0.81 

16.02 

0.160 

59.1 

1.22 

2.91 

0.041 

10.7 

9.47 

28.43 

0.262 

105.0 

0.63 

26.52 

0.240 

98.1 

0.67 

21.64 

0.198 

79.9 

0.83 

11.54 

0.144 

42.5 

2.11 

7.39 

0.099 

27.1 

3.55 

Table  3  (cont’d), 


Cix/0  ‘ 
tg/cm  '  sj 

f-Ap ), 
Igf/eml 

Re 

fp 

/2g,tV-Ap7,i'/r(ic;7 

5.02 

0.074 

18.4 

5.73 

9.52 

0.127 

35.0 

2.74 

29.18 

0.319 

107.8 

0.73 

22.57 

0.243 

83.3 

0.93 

19.32 

0.211 

71.2 

III 

34.16 

0.327 

126.3 

0.55 

30.59 

0.290 

113.0 

0.61 

27.75 

0.265 

102.5 

0.67 

14.43 

0.149 

53.2 

1.40 

15.43 

0.150 

56.9 

1.23 

25.50 

0.212 

94.2 

0.64 

26.98 

0.226 

99.7 

0.61 

24.68 

0.202 

91.2 

0.65 

6.90 

0.096 

25.4 

3.94 

22.68 

0.321 

83.8 

1.22 

9.19 

0.135 

33.8 

3.12 

11.91 

0.172 

43.9 

2.37 

7.98 

0.121 

29.4 

3.71 

19.77 

0.266 

73.0 

1.33 

21.74 

0.290 

80.3 

1.20 

27.53 

0.357 

101.8 

0.92 

22.46 

0.271 

83.0 

1.05 

23.41 

0.283 

86.5 

1.01 

26.10 

0.310 

96.6 

0.89 

27.37 

0.320 

101.3 

0.84 

28.64 

0.311 

106.0 

0.74 

29.31 

0.316 

108.5 

0.72 

32.42 

0.313 

120.1 

0.58 

36.38 

0.340 

134.8 

0.50 

38.82 

0.358 

144.0 

0.46 

26.10 

0.221 

96.6 

0.63 

20.61 

0.219 

76.2 

1.01 

3.02 

0.035 

III 

7.50 

18.97 

0.202 

70.1 

1.10 

20.60 

0.220 

76.2 

1 .01 

23.36 

0.228 

86.4 

0.82 

17.68 

0.164 

65.3 

1.03 

25.68 

0.236 

95.1 

0.70 

26.73 

0.248 

99.0 

0.68 

28.09 

0.238 

104.0 

0.59 

28.82 

0  250 

106.7 

0.59 

30.94 

0.258 

114.5 

0.53 

24.60 

0.201 

91.0 

0.65 

G  x  /0  *  f-Ap ),  Re  fp 

(g/vm1-.)  (gf/cm!J  rD„FfeG/,.y  /2g,  D^-Ap/^/F^G'/ 


e«  =  0.447  g/em‘,  l)r  =  0.092  cm 


2.38 

0.053 

6.2 

13.05 

3.59 

0.089 

9.3 

9.62 

8.00 

0.166 

20.8 

3  60 

32.03 

0.599 

83.8 

0.81 

30.71 

0.566 

80.3 

0.83 

29.87 

0.544 

78.2 

0.84 

28.64 

0.507 

75.0 

0.86 

24.82 

0.425 

64.9 

0.96 

13.34 

0.224 

34.8 

1.75 

16.93 

0.281 

44.2 

1.37 

qs  -  0.472  g/em 

\  l)p  =  0.065  cm 

7.79 

0.379 

14.3 

6.13 

13.68 

0.667 

25.2 

3.50 

24.30 

1.120 

44.8 

1.86 

13.70 

0.617 

25.2 

3.23 

22.31 

0.909 

41.1 

1.79 

14.60 

0.553 

26.9 

2.55 

2.84 

0.140 

5.2 

17.06 

1.62 

0.098 

3.0 

36.89 

21.78 

1.142 

40.2 

2.36 

34.31 

1.755 

63.5 

1.46 

40.89 

2.122 

75.7 

1.25 

40.75 

2.016 

75.4 

1.19 

38.14 

1.871 

70.6 

1.26 

10.41 

0.526 

19,1 

4.77 

44.89 

2.098 

83.1 

1.02 

42.66 

1.976 

79.0 

1.07 

40.72 

1.868 

75.4 

111 

7.65 

0.364 

14.1 

6.10 

45.72 

2.020 

84.7 

0.95 

30.26 

1.320 

56.0 

1.42 

31.10 

1.326 

57.5 

1.35 

31.58 

1.348 

58.4 

1.33 

67.45 

2.862 

125.2 

0.62 

52.32 

2.175 

97.1 

0.78 

17.90 

0.710 

33.0 

2.18 

47.02 

1.828 

87.1 

0.81 

45.59 

1.753 

84.5 

0.83 

49.67 

1.879 

92.0 

0.75 

55.40 

2.066 

102.8 

0.66 

fj  "(q/cm3) 

<•)  0.377 
M  0.387 
(0  0.427 
(•)  0.447 
(*)  0  472 
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Figure  3.  Relationship  between  friction  factor  fp 
and  Reynolds  number  Rep. 
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Figure  4.  Relationship  between  friction  factor  Vp 
and  Reynolds  number  ReL 
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As  shown  in  Table  2,  the  pressure  drop  through 
the  snow  layer  for  snow  with  density  ranging  from 
0.377  to  0.472  g/cm1  can  be  fairly  well  represented 
linearly  with  the  superficial  mass  flow  rate  G.  The 
least-squares  analysis  for  each  snow  density  result¬ 
ed  in  correlation  coefficients  from  a  moderate  val¬ 
ue  of  0.8130  to  a  rather  high  value  of  0.9933. 
From  Figure  3,  it  can  be  seen  that  the  air  flow  en¬ 
counters  higher  resistance  through  snow  as  com¬ 
pared  to  gas  flow  through  randomly  packed  por¬ 
ous  columns  of  metallic  and  nonmetallic  spherical 
and  nonspherical  materials  (in  the  latter  case, 
there  was  no  phase  change  or  mass  transport  ac¬ 
companying  the  flow).  Since  snow  is  very  sensitive 
to  temperature,  its  physical  properties  are  con¬ 
stantly  changing  through  metamorphism,  which 
creates  ice  bonds  among  the  snow  grains  even 
under  isothermal  conditions  through  the  process 
of  sintering. 

Since  the  experiment  was  set  up  to  evaluate  the 
thermal  conductivity  and  vapor  diffusivity  of  the 
snow,  a  moderate  temperature  gradient  was  im¬ 
posed,  even  though  the  air  entering  at  the  top  of 
the  column  was  saturated  to  prevent  or  to  limit  the 
mass  transport  between  the  snow  grains  and  the 
fluid  stream  as  it  flowed  through  the  interconnect¬ 
ed  pores.  Though  a  thermal  equilibrium  is  as¬ 
sumed  between  the  snow  grain  and  the  air  stream 
as  it  passes  the  grain  surface,  the  air  still  picks  up 
the  water  vapor  as  it  flows  downward  and  comes 
into  contact  with  snow,  which  was  maintained  at 
higher  temperature  at  the  exit.  The  moist  snow 
must  exhibit  more  drag  or  stickiness  to  the  air  flow 
than  the  dry  and  rather  smooth  surfaces  of 
metallic  and  nonmetallic  surfaces. 

In  laminar  flow,  equations  of  the  same  form 
may  be  used  for  the  friction  loss  through  the  por¬ 
ous  medium  as  for  the  fluid  flow  through  the  con¬ 
duits  (i.e.  the  pressure  drop  due  to  friction  losses 
through  a  porous  medium)  is 


/p  - 


Rep 


The  only  difference  between  /p,  Rep  and  /,  Re  is 
that  /p  and  Rep  contain  factors  of  F,  and  Fke,  re¬ 
spectively,  in  addition  to  the  parameters  used  in 
defining/and  Re.  By  rewriting  eq  21,  i.e.  by  solv¬ 
ing  for  velocity  u,  it  follows  that 


Sc^p^Rc  Q<w  f  K  Q>w  f 

32  F,  Ln  Ln 


(24) 


in  which  K  is  the  air  permeability  of  snow  used  in 
eq  1,  i.e. 


K  —  K< 

K  32  Ff 


(25) 


From  this,  the  value  of  K  can  be  calculated  with¬ 
out  performing  an  experiment,  if  the  dimension¬ 
less  factors  ^Rc  and  Ff  and  the  particle  dimensions 
are  known.  The  value  of  K  is  directly  proportional 
to  Dp  but  implicitly  K  depends  on  porosity  and 
sphericity  through  the  factors  of  ^Re  and  Ff. 

Bender  (1957)  reported  an  extensive  study  on  air 
permeability  of  snow  and  pointed  out  that  the  air 
flow  is  laminar  for  velocities  of  less  than  5  cm/s 
for  fine-grained  snow  (less  than  0.8  mm  in  diam¬ 
meter);  2  cm/s  for  medium-grained  snow  (0.8  to 
1 .2  mm  in  diameter);  and  1  cm/s  for  large-grained 
snow.  He  reported  that  all  of  his  results,  as  long  as 
the  flow  is  in  laminar  region,  can  be  represented 
by  a  single  relationship,  i.e. 


K  = 


16.80^,  (0 

fo  -  f 


(26) 


where  K  =  air  permeability  expressed  in  cm/s 
Dp  -  average  snow  grain  size  in  mm 
f0  =  highest  porosity  (i.e.  when  the  packing 
is  in  loosest  state)  attainable  for  the 
given  shape  of  the  snow  particle. 


(-Ap)f  = 


ttLvuFf 

gcDj,FR' 


e/Wf 


and  for  flow  through  a  conduit  is 


(21) 


(-Ap)f  =  =  elw,.  (22) 

In  both  flows,  the  friction  factor  can  be  represent¬ 
ed  by 


Since  no  specific  experimental  data  were  given  in 
terms  of  pressure  gradient,  volumetric  flow  rate  or 
superficial  ’  elocity,  or  porosities  t  and  r0,  it  was 
not  possible  to  compare  the  results  from  eq  25  and 
26.  Furthermore,  snow  is  so  sensitive  to  the  tem¬ 
perature  and  aging  through  metamorphism  that  it 
will  be  indeed  difficult  to  reproduce  experimental 
results  even  if  identical  samples  and  measurement 
techniques  are  used. 
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It  can  be  concluded  from  this  study  that  the 
pressure  drop  or  friction  loss  of  air  flow  through  a 
uniform  snow  layer  can  be  correlated  in  terms  of 
/p  and  Rep  with  a  high  correlation  coefficient 
(0.9688)  or  in  modified  form  of  fj  and  Rep  with  a 
moderate  correlation  coefficient  of  0.8385.  In 
view  of  the  complex  processes  of  mass  exchange 
and  physical  structure  transformation  taking  place 
even  under  isothermal  conditions,  the  results  pre¬ 
sented  should  provide  useful  information  in  pre¬ 
dicting  energy  requirements  for  air  to  penetrate  a 
uniform  snow  layer  of  various  snow  densities. 
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